Braided Bundles and Compact Coils: The Structure and Thermodynamics of 
Hexagonally-Packed, Chiral Filament Assemblies 
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Molecular chirality frustrates the two-dimensional assembly of filamentous molecules, a fact that 
reflects the generic impossibility of imposing a global twisting of layered materials. We explore the 
consequences of this frustration for hexagonally-ordered assemblies of chiral filaments that are finite 
in lateral dimension. Specifically, we employ a continuum-elastic description of cylindrical bundles 
of filaments, allowing us to consider the most general resistance to and preference for chiral ordering 
of the assembly. We explore two distinct mechanisms by which chirality at the molecular scale of the 
filament frustrates the assembly into aggregates. In the first, chiral interactions between filaments 
impart an overall twisting of filaments around the central axis of the bundle. In the second, we 
consider filaments that are inherently helical in structure, imparting a writhing geometry to the 
central axis. For both mechanisms, we find that a thermodynamically-stable state of dispersed 
bundles of finite width appears close to, but below, the point of bulk filament condensation. The 
range of thermodynamic stability of dispersed bundles is sensitive only to the elastic cost and 
preference for chiral filament packing. The self-limited assembly of chiral filaments has particular 
implications for a large class of biological molecules - DNA, filamentous proteins, viruses, bacterial 
flagella - which are universally chiral and are observed to form compact bundles under a broad 
range of conditions. 

PACS numbers: 87.16.Ka, 61.30.Dk, 87.15.bk 



I. INTRODUCTION 



From a sufficiently close perspective, understanding 
the effect of chirality on the organization high-aspect 
ratio molecules appears to be problem of local geom- 
etry. The lack of mirror symmetry generically implies 
that rod-like molecules exert a mutual torque, favoring 
a preferred tilt of neighboring molecules as opposed to a 
parallel arrangement of molecular axes favored in achiral 
systems It is this vantage point that originally com- 
pelled Crick to propose the so-called "coiled-coil" model 
to explain the structure of a-karetin from x-ray diffrac- 
tion data But it is well-known that chiral forces, 
while local in nature, have profound implications for the 
global arrangement of molecules in ordered systems. A 
fascinating example of this global influence occurs in the 
"blue phases" of liquid-crystals, wherein mesoscopic chi- 
ral order is delicately maintained by a periodic network 
of topological defects Q. A more explicit example of 
this effect was noted by de Gennes for chiral liquid crys- 
tals with smectic order Q. Just as the application of 
magnetic field destroys the superconducting state, the 
tendency for chiral, or cholesteric, ordering disrupts the 
ability to form a periodic stack of layers. At best, chiral 
order exists in a smectic system in the neighborhood of 
screw-dislocation grain boundaries, where periodic order 
has locally been destroyed [f|. 

Beyond the one-dimensional case of smectic materials, 
chiral-tilt order is more generally incompatible with bulk 
periodic order. Indeed, it routinely observed that bulk 
systems of chiral polymers-such as DNA fragments-expel 
twist above a critical concentration at which they adopt 



two-dimensional, columnar, liquid-crystalline order [y]. 
For hexagonally-ordered chiral polymers, tilt order that 
reflects that chirality is only possible in the presence dis- 
location grain boundaries that locally tear the underly- 
ing two-dimensional lattice 0, [H ■ That columnar phases 
have this property in common with smectic phases is not 
a surprise, as the array of two-dimensionally ordered lines 
can be decomposed into layers, each of which is incom- 
patible with a global twist. 

While global chiral ordering is not possible for a bulk 
periodic system, the same may not be said for a finite do- 
main. This fact is well-appreciated in the context chiral 
membranes 0, [ToL [Tlj , since ribbon of finite width may 
be twisted into a helicoid without tearing. Hence, the 
material is globally ordered in a chiral sense, while at the 
same time maintaining the inherent layer geometry. In 
a recent Letter [12| , we reported on a similar mechanism 
in which molecular chirality induces a global twist of a 
two-dimensional assembly of filaments, provided that the 
assembly is finite in diameter. The chiral ordering of fila- 
mentous assemblies has profound implications, with par- 
ticular consequences for systems of biological filaments, 
such as DNA or filamentous proteins. Namely, global 
twisting of bundles can ultimately limit the stable size 
bundle aggregates, providing a generic and robust means 
by which cells may regulate the assembly of a host of 
biological filaments into fiber bundles. It is the purpose 
of this article to present and analyze a general model of 
filament assembly in which molecular chirality frustrates 
the ability to form two-dimensionally ordered aggregates, 
ultimately leading to a state of thermodynamically stable 
bundle assemblies of finite diameter. 

In this article, we explore two geometrically distinct 
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FIG. 1: Two packing motifs through which molecular struc- 
ture of constituent filaments imparts chirality on the orga- 
nization of densely packed bundles. In (a), intermolecular 
forces of chiral filaments favor relative tilting of neighboring 
filaments, leading to an overall twisted bundle. In (b), the 
unstressed shape of the filament itself is helical, imparting a 
writhe to the central axis of the bundle. 



mechanisms through which chirality at the molecular 
scale frustrates the two-dimensional assembly of fila- 
ments. In the first mechanism, chirality is a colligative 
effect, favoring a difference in backbone orientation when 
two filaments are brought into close contact. Locally, fil- 
aments adopt the coiled-coil packing preferred by inter- 
locking molecular screws. As a filaments are added to a 
growing bundle, they twist, or braid, helically around the 
central axis of the bundle (see Figure [2(a)), adopting the 
so-called "double-twist" geometry of blue phase liquid 
crystals fl3| . Such a mechanism was originally proposed 
to explain the observed twisting of fibrin bundles [lj| , as 
well as to explain two-dimensionally ordered and twisted 
structures observed in chromatin of dinoflagellates 
In ref. fl2l ]. we established the thermodynamic viabil- 
ity of this mechanism for the first time and explored the 
sensitivity of the limited bundled growth to exact type 
and quality of two- and three-dimensional order in the 
assembly. 

Here, we present a more detailed analysis of the model 
of twisted, straight bundles. We find that the optimal 
degree of bundle twist is a generically non-mononotonic 
function of bundle radius, with a single maximum whose 
location is determined only by the mechanical resistance 
of the bundle to twist. In the simplest case, that of a 
hexagonal-columnar, liquid-crystalline bundle, this char- 
acteristic size is the bend penetration depth, cx 
(K 3 / fi^) 1 ^ 2 , where K s and are the respective elas- 
tic constants describing the resistance to bending and in- 
plane shear distortions of the hexagonal filament packing. 



The predicted non-monotonic dependence optimal twist 
on bundle radius is critically linked to appearance of an 
optimal bundle radius, which is in the range of the char- 
acteristically mesoscopic length scale, A3^. We find that 
a dispersed state of finite-sized bundles is predicted to be 
stable for systems near to, but below, the point of bulk 
condensation, in which aggregate size is thermodynami- 
cally unlimited. 

In addition to the case of twisted bundle, here we ex- 
plore a novel mechanism in which molecular chiral struc- 
ture obstructs the ability to form two-dimensional fil- 
ament assemblies. In this second mechanism, chirality 
enters through the inherent helical structure of filaments 
themselves. This is the well-known geometry of bacterial 
flagella, whose "cork-screw" morphology sensitive to so- 
lution conditions as well as to mechanical stress [ijj • 
There is also some evidence to suggest that certain rod- 
like viruses, like the fd virus, adopt a weakly supercoiled 
geometry in presence of certain modifications of the pro- 
tein coat fig . [HI . It is already clear from observations 
of liquid-crystalline mesophases of the salmonella flagella 
that helical backbone structure has a dramatic influence 
on the intermolecular packing of filaments [2(j. In the 
two-dimensional filament packings considered here, the 
structure of the individual filament imparts an overall 
writhing geometry to the central axis of the bundle (see 
Fig. [2(b))- We consider the case of isometric packings of 
filament recently analyzed in detail by Starostin, in which 
the bundle maintains a perfectly hexagonal arrangement 
perpendicular the filament backbone [211] . While these 
bundles are optimal in terms of constant preferred dis- 
tance between neighboring filaments, growing a bundle to 
finite radius requires filaments at the periphery of the to 
be distorted from their preferred state of writhe. Thus, 
we find that the preferred geometry of filaments also frus- 
trates the growth of hexagonally-packed bundles of heli- 
cal filaments. 

Although this second mechanism is quite distinction 
from the case where filaments twist around a straight 
bundle axis - isometric bundles have strictly no twist - 
the thermodynamic consequences are quite similar. For 
solutions of helical filaments sufficiently close to, but be- 
low, the point of bulk condensation, bundles of filaments 
grow to finite diameter. As the bundle grows, the bun- 
dle unwinds, straightening out the helical filaments. As 
such, we find that the range of thermodynamic stabil- 
ity of finite-sized helical bundles is sensitive only to the 
mechanical cost of unbending the filament from its pre- 
ferred state. That is, in order to form a bulk aggregate, 
the net cohesive energy gain per bundled filament, must 
be larger than the mechanical cost of straightening a he- 
lical filament. 

This paper is organized as follows. In Sec. [IT] we 
present a generic, continuum elastic model for hexago- 
nal assemblies of chiral filaments. In Sec. IIIII we analyze 
the thermodynamic behavior of hexagonal bundles in the 
presence of twist-inducing interactions between chiral fil- 
aments. In Sec. IIVI we analyze a model of helical fila- 
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merits that form hexagonally-packed, writhing bundles. 
We conclude with a brief discussion in Sec. [V] 



II. ELASTICITY OF HEXAGONALLY-PACKED 
FILAMENTS 

The structure and mechanical properties of a filament 
bundle derive from the microscopic interactions between 
filaments. While these microscopic interactions may be 
quite complex, we may understand the inherent cost 
of long-wavelength distortions of the bundle by consid- 
ering generic form of the elastic Hamiltonian for two- 
dimensionally ordered filaments. These most general de- 
scription of the elastic properties can be written in terms 
of u_i_ (x) , which describes the local disp lacement of a fila- 
ment from its equilibrium position [13| . For hexagonally- 
ordered filaments aligned along the z axis, the elastic en- 
ergy can be constructed from the in-plane strain tensor, 
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(diU±j + dju± t - d l u± ■ diU± - Utj) . (1) 



Here, i and j refer only to in-plane directions (aij _L z) 
and t ~ z + d z u± is the unit tangent vector describing 
filament orientation. The first non-linear term in uh is 
the same non-linear correction appearing for purely two- 
dimensional solids, while the second non-linear term, fa- 
miliar from the non-linear elastic theory of smectic liquid 
crystals, ensures rotational invariance of the elastic en- 
ergy around in-plane axes [13] • We will see below that it 
is necessary to retain the non-linear energy contributions 
resulting from these "gauge" terms. 

In the absence of broken translational symmetry along 
the filament axis, the array responds to elastic deforma- 
tions as a columnar-hexagonal liquid crystal, 



(2) 



where A_l and fi± are the Lame constants describing 
the resistance to compressive and shear distortions of 
the hexagonal order perpendicular to the filament axis. 
Higher-order derivative contributions take the form the 
of the Frank free energy for polymer nematics, 

Hi = iyd 3 .T{Xi(V ± -t) 2 +^ 2 (t-Vxt) 2 +X 3 |(t-V)t| 2 

+ K 2i V- [(t- V)t-t(V-t)]}, (3) 

where K\ , K 2 and K3 are the respective elastic constants 
to splay, twist and bend of the filament backbone. The 
total-derivative term on the second line of ((3]) is the so- 
called "saddle-splay" term, which is often neglected for 
bulk systems [13[ . For a typical 3D solid, one may ignore 
the higher-order Frank free energy contributions of eq. 
(J3j> , but in the absence broken-translational symmetry 
along the backbone direction it is necessary to include 
at least the resistance to bending (K3 7^ 0) in order to 



stabilize the system to long- wavelength deformations and 
thermal fluctuations. 

Below a critical temperature, a columnar system of 
filaments will pass to a crystalline state [H|,[24j and adopt 
the elastic response of a 3D anisotropic solid. The strain 
energy of a solid is properly described by the 3D strain 
tensor 25], 



-(diUj 



djUi 



diU ■ dju), 



(4) 



where u here is the three-dimensional displacement and 
i and j refer to all three spatial coordinates. The most 
important distinction between 3D-solid and 2D-columnar 
response is a resistance to uniform shear along the fila- 
ment axis. In addition to the elasticity described by eq. 
@ a 3D-solid requires the additional elastic terms, 
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d 3 x{\\ l u 2 zz + 2n ll (u 2 xz + u 2 yz )}, (5) 



where, to be clear, this contribution to the elastic re- 
sponse is written in terms of the 3D strain tensor, not 

u 1 -- 
u v . 

The elastic response described by Tt±_ + Tin, provides 
a fully general description of the a 3D hexagonal solid. 
By taking either 7^ or =0 we may use this 
description to model a bundle with either solid or liquid- 
crystalline order. Alternatively, we may view the com- 
posite elastic response in terms the following microscopic 
model. Consider a hexagonal array of inhomogeneous 
filaments where intermolecular forces depend indepen- 
dently on (1) the in-plane separation between neigh- 
boring filament backbones and (2) the separation be- 
tween mass points distributed inhomogeneously along 
the filament backbone. The former responses derives, 
say, from the monopole interaction between charged fil- 
aments, while the second results from the combination 
of multi-pole electrostatic, steric and protein- or ion- 
mediated forces at work in biofilament systems. For 
hexagonally-ordered filaments we must always have Aj^ =/= 
and ^ 0, while for sufficiently low concentrations 
of linking agents or sufficiently strong thermal fluctua- 
tions /iy may be considerably reduced relative to the in- 
plane moduli. In the following, we treat these moduli 
as mesoscopic parameters which may be tuned through 
underlying microscopic physics. 

Finally, we consider the elastic terms that reflect 
the microscopic chiral structure of biological filaments. 
These terms break the manifest chiral symmetry of eqs. 

^ and ([5]) and change sign under spatial inversion, 
x — * — x, but do not change sign under t — > — t. For fila- 
ments possessing hexagonal order the three most relevant 
chiral elastic terms are, 

H* = J d 3 x{ji • (V X t) + 7'(t • V) [t • (V x u x )] 

+ 7 "t- [(t-V)Vxux]}. (6) 



When uj^ is small, each of these terms is linear in uj_ 
and second order in derivatives. The first term is simply 
the cholesteric twist of the filament backbone typically 
appearing for chiral nematic liquid crystals [l3| . The sec- 
ond and third terms are rotationally-invariant generaliza- 
tions of the a chiral symmetry-breaking term appearing 
in columnar phases, studied by Kamien and Nelson in 
the context of chiral defect phases [7], [8| . To first order 
in u^, each of these terms are equivalent and related to 
the rotation of the hexagonal bond orientation along the 
column axis, 8 z 0q, where 9q is the six- fold bond angle 
of a hexagonal lattice measured with respect to a refer- 
ence direction in the plane. Hence, for filament bundles, 
chirality of the filaments induces "braiding" of filaments 
around the bundle axis. Pursuing these terms to higher 
order, we note that the final two terms in (JSJ) are different 
since 

V x u ± = 28 6 i + 1 x (t • V)u ± 

~ (V x ux) z z + z x d z uj_, (7) 

The second term in ^ only picks up <9 Z (V x ujj 2 while 
the third has an additional contribution, d z u± ■ (z x 
d z u±). This higher order contribution is a measure of the 
writhe of the filament backbone, representing an impor- 
tant class of geometrically distinct chiral deformations. 
Specifically, it may be shown that to lowest order in , 
all three terms in eq. ([6]) are proportional to the local 
twist of filaments around the central axis of the bundle, 
while d z u± ■ (z x d z u±) is proportional (to lowest or- 
der) to the writhe of the backbone of the bundle. Twist 
and writhe terms are more commonly used to describe 
the geometry of linked curves, such as the two nulceotide 
strands of a DNA molecule [49|. In the Appendix, we 
show that relationship between the chiral elastic terms 
of a hexagonally-ordered filament array in eq. ([6]) can be 
more formally identified with either the bundle twist, or 
the local contribution to filament writhe. 

In the following sections, we treat these two effects 
separately. In Sec. IIII[ we consider the effect of twist- 
inducing terms, like t • (V x t), in hexagonal bundles of 
otherwise straight filaments. There is an important dis- 
tinction between the twist- and writhe- inducing terms. 
The twisting of filament measures differences in geome- 
try (orientation) between distinct filaments in the bun- 
dle. Therefore, it is a colligative affect, reflecting the 
chiral nature filament interactions. On the other hand, 
the writhe is really a measure of differences in geometry 
along a single filament (see eq. IA.6|) . The latter class of 
energetic terms would be present even in the absence of 
neighboring filaments. Unlike the twist-inducing terms, 
the writhe-inducing terms in H.* is really a measure of the 
intrinsically preferred geometry of the filaments them- 
selves. Therefore, we take an alternative approach to 
study the assembly of bundles induced to writhe. In Sec. 
IIVI we consider hexagonal assemblies which are induced 
to writhe by preferred helical structure of the filaments 
themselves. 
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FIG. 2: (a) shows a cartoon of the torsional deformation of a 
vertical section of the bundle near z = 0, demonstrating that 
the hexagonal packing is largely preserved. The filled (open) 
circles show the original (deformed) position of filaments, (b) 
depicts cylindrical sections from a twisted bundle at various 
radii. Starting from the center and moving outward, accord- 
ing to eq. (|I0p . the tilt of filaments relative to the central axis 
grows approximately linearly. The tilt angle of outermost fil- 
aments is then 6 = QR, which leads to a reduction of the 
spacing between adjacent filaments, d, below the preferred 
spacing, do- 



III. STRAIGHT BUNDLES OF TWISTED 
FILAMENTS 

In this section, we analyze the elastic energy intro- 
duced in Sec. HH for bundles whose response to the in- 
trinsic chiral stress eq. ^ is to braid around a central 
axis that is uniformly straight. This case that was dis- 
cussed in ref. [12], in the context of self-limited growth 
of biofilament bundles. The torsional deformation of a 
cylindrical bundle is described by the following in-plain 
displacement (in cylindrical coordinates) 

u ±( x ) = r[cos(ilz) — llf + r sm(Qz)<f), (8) 

where 2tt/H, corresponds to the pitch of filaments winding 
helically around the z. axis. From eq. (|TJ) it is straightfor- 
ward to show that this distortion leads to the following 
in-plain strain, 

hence, the strain-energy grows radially as (fir) 4 . Geo- 
metrically, this strain energy results from the rotation of 
the filament tangent, relative to the z axis, 

i~z + rir$. (10) 
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Figure [2] shows a bundle whose outer filaments have been 
rotated to an angle 9 ~ ilR relative the central axis of 
the bundle. The distance between filaments along the 
azimuthal direction is the same as preferred filament sep- 
aration in the center of the bundle, cIq. Thus, separation 
of filament perpendicular to the filament axis at the outer 
edge of bundles is reduced according tod 2 ~ d 2 , ( 1 — 9 2 /2 ) , 
which indicates that compressive and shear strains of or- 
der 9 2 are introduced by the torsional strain. 

Depending on the relative magnitude of the in-plane 
Lame constants, the elastic energy may be lowered by 
superimposing the displacement of eq. (jS| with an ad- 
ditional dialation in the radial direction according to 
U J-( X ) ~~ * U J-( X ) + (o:n 2 r 3 /2)f, where a is variational 
parameter. This deformation gives the following contri- 
bution to the in-plane elastic energy (correct to 0[(f2r) 6 ] 

), 



Ha 



+ ^(l + 2a) 2 ]}. (11) 



The variational parameter, a, is determined by minimiz- 
ing the elastic energy and hence this parameter depends 
only on the relative magnitude of the Lame coefficients. 
In the limit that the bundle is incompressible in the plane 
and A_l 3> fJ-±, a = 1/3 and for this torsional defor- 
mation the elastic due hexagonal order of filaments is 
H±(Q)/V = (25/216)/i ± (fyT) 4 , where V is the volume 
of the cylindrical bundle [51|. Note that no choice of 
a allows us to relieve the elastic stress introduced by 
bundle twist. The elastic strain introduced by a double- 
twist configuration of lines is a geometric consequence of 
non-zero saddle splay (26j . which is incompatible with a 
constant 6-fold coordinated lattice geometry. 

For filament bundles with solid elastic response, the 
torsional displacement of eq. ([8]) also leads to shear con- 
tributions to the out-of-plane and divergence-free stress 
contributions, 



a xz = /i||Oy; a yz 



-/i. 



Qx. 



(12) 



The radial force due to these stresses on the boundary of 
a cylindrical bundle vanishes, as required by mechanical 
equilibrium. The linear-elastic response of a solid rod 
to a torsional deformation is well known, Wii(f2)/V = 
/i||(flR) 2 /4 [|3. 

The Frank-energy contributions can be computed di- 
rectly from eqs. and (JSJ) to show, 



Mt(") 
V 



= K 2 - 



K 2i 

2 



n 2 



Or 



(13) 



Here, we note that the sign of K 24 is not constrained by 
symmetry so that for sufficiently large saddle-splay con- 
stants (K 24 > 2K 2 ), the Frank elastic energy is unstable 
to twist, even for achiral systems. Indeed, model calcu- 
lations for achiral carbon nanotube ropes suggest that 



highly adhesive van der Waals interactions drive sponta- 
neous braiding of nanotube ropes [28| . Finally, we sum- 
marize the contribution from the chiral-elastic terms in 
the bundle as H*(Q)/V = -2-f Tw n(l + 0[(nR) 2 }) , where 



7Tw = 7 + l' + i ': 



(14) 



is the total coupling to the local twisting of filaments 
around the central axis of the bundle (see Appendix) . 

The total elastic energy of the twisted bundle is the 
sum of these components discussed above, E tW i S t{P) = 
H±(fl) + H\\(n) + HiiCl) + H*(Q). By defining the fol- 
lowing length scales, 



\2 _44X 3 . 
As± = 25^? 2 ' 



^ 2(2K 2 - K 24 ) 

Mil 



and angles 



2 _ 2(2K 2 - K 2i ) 
?23 - 77 ; 



_ 44 M || 
~ 25^i± 



(15) 



(16) 



the clastic energy for a straight bundle of chiral filaments 
takes on the general form, 



Ei 



. R) 



nKvL 



/ R 2 



R 2 \ 
^)-l9R 

where 9 = CIR is the twist angle of outermost filaments 
in the bundle (see Fig. [2]) and 7 = 87t w /^3- For given 
bundle radius, the equilibrium twist angle derives from 
the solution to dE tW i S t/d9 — 0. This solution to this cu- 
bic equation predicts the full dependence of bundle twist 
on the many elastic parameters in model. Many of these 
phenomenological elastic parameters may be difficult ei- 
ther to predict from microscopic considerations or to de- 
termine independently from macroscopic measurements 
of filament assemblies. We focus here on the dependence 
equilibrium bundle twist, 9, on the radius of the bundle, 
which may be tested directly by structural observations 
of bundles under various conditions. 

Despite the many parameters that specify the elastic 
energy, the mechanical equilibrium allows for a rather 
limited range of radial dependence of twist. The chiral 
term always drives # ^ 0, while the remaining terms re- 
sist torsion. A simple inspection of eq. (TIT)) demonstrates 
that the form of the dominant elastic restoring force, be it 
twist, bend or shear stiffness, depend only on the values 
of 9 and R relative to the angle and length scales defined 
in eqs. (fT5"|) and (flB)) . The full dependence minima of 
Etwist on R is shown in Figure [3) 

For small bundles, twist generically grows with bun- 
dle radius: for large (small) torsion, bundle twist is 
restrained by the Frank bend (twist) response. While 
for large bundles, equilibrium twist generically decreases 
with bundle size: at large (small) 9, bundle twist is lim- 
ited by out-of-plane (in-plane) resistance to shear. As 
depicted in Fig. [3J the growth of 9 for small bundles and 
subsequent decrease for large R only allows for a single 
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e = Sin 



FIG. 3: Curves predicting the radial-dependence equilibrium 
twist angle for fixed elastic constants-according to the min- 
ima of eq. (|17p ~are shown in (a) and (b). Each solid curve 
depicts a different value of the preference from twist 7tw 
Depending on the relative sizes of elastic parameters A3X and 
A211, the variation of 6 with R may be summarized by one of 
two scenarios, which are shown in (a) and (b) and described 
by 4 possible scaling regimes. Note, in particular, that 9{R) 
has only a single maximum. 



maximum in the curve 0{R). The non-monotonic de- 
pendence of induced twist on bundle size can be directly 
attributed to the self-limited thermodynamic growth of 
chiral filament bundles. The thermodynamic behavior of 
this model has been summarized in a previous report [l2| . 
We therefore present here a brief discussion of two cases: 
assembly of columnar-liquid crystalline bundles and as- 
sembly of strongly cross-linked, solid bundles. 



A. Columnar Liquid-Crystalline Bundles (piy = 0) 

Above a critical temperature or below a critical den- 
sity, filaments in a hexagonal bundles are free to slide lon- 
gitudinally with respect to their neighboring filaments. 
This corresponds to the columnar liquid-crystalline state 
of filament order characterized by a vanishing resistance 
to shears along the filament axis =0). In the model 
discussed above, this corresponds to the case of = 0. 
For the sake of simplicity, we also assume in the follow- 



ing that 6*23 = 0. The physical basis of this assumption is 
that the bend elastic constant K3 is determined primar- 
ily by the intrinsic stiffness of the filaments themselves, 
which we expect to be sufficiently large that K2/K3 <C 1. 

In this case the elastic energy for a twist bundle takes 
the form, 



E,, 



I R) K; 



3M 



irL 



R 2 

A 3± 



2 lTw 6R, (18) 



For this energy it is straightforward to compute the pre- 
ferred value of bundle twist, Oq, for a fixed size, 



2rrwV/3 



R 1 / 3 



(l + Ryxl^/f 



(19) 



Note, specifically, the non-monotonic dependence of 
8q(R). There is a single maximum at R = A3j_, and 
6 (R) ~ i? 1 / 3 for R < A 3 _l while 6 (R) ~ iT 1/3 for 
R ^3±- Physically, this crossover is due to the fact 
that bending energy dominates the mechanical resistance 
to twist for R < A3 , while the in-plane shear dominates 
in the opposite regime. Using this result in E twist (9o, R) 
we may compute the radial dependence of the energy 
gained by the induced twisting of the bundle, 



E, 



twist ("0, 



R) 



Tw 



7rL 



3/2^ 
2 V K 3 



1/3 



i? 4 / 3 



(1 



3± 



)l/3- 



(20) 



Since E tW i S t(&o, R) ~ —OR, this means that the elastic 
twist energy gain grows faster than linearly, — i? 4 / 3 , for 
small bundles and more slowly than linearly, — i? 2 / 3 , for 
large bundles. This result means that the lateral force 
on the bundle due to elastic forces, —dE tW i St /dR, van- 
ishes in both the small bundle and large bundle limits. 
Under appropriate thermodynamic conditions, this prop- 
erty gives rise to a thermodynamically preferred size, Rq, 
which is finite. 

To see this, we consider additional contributions to the 
free energy of a growing bundle that are due to cohesive 
forces between neighboring filaments. In general, we ex- 
pect that these forces can contribute a net free energy 
gain per bundled filament, — AeL. To be clear, Ae is 
the cohesive energy gain per unit length of parallel fila- 
ments in a hexagonal bundle. Due to the chiral nature 
of filament interactions, the twisting of the bundle may 
further increase the cohesive energy of the bundle, and 
this mechanism is described by E tW i S t ■ Additionally, due 
to the reduced number of adhesive contacts and the edge 
of the bundle, we may also attribute a positive energy 
per unit area of the bundle surface, S. 

Combining these twist-independent contributions with 
Etwist we have the total free energy of a bundle, 



F(R) 



= 2irER - nR 2 p Ae 

37r/2 7 4 w y/3 i? 4 / 3 

2\K 3 ) (1 



3± 



)V3 



(21) 
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where po is the number of filaments per unit area of the 
bundle cross-section. 

At equilibrium, aggregates of radius R exist 
in solution according to a probability p{R) oc 
exp{— F(i?)/fcsT} [30j. The dominant state of the so- 
lution is determined by the global free energy minimum 
of F(R). Respective minima at R — > and R — > oo 
correspond to the states of dispersed single filaments in 
solution and bulk aggregates of filaments of macroscopic 
size. In the absence of any twist-dependent energy gain 
(7Tw — > 0) , a state of bundles of finite radii dispersed in 
solution does not occur. 

When Ae > 0, then the bundle free energy is un- 
bounded at R — ► oo, indicating a state of bulk aggre- 
gation (infinite aggregates). However, when the net co- 
hesive energy for the parallel state of neighboring fila- 
ments is non-positive, two possible states arise. This is 
most clear from the case when Ae = in F(R). The 
positive surface energy grows linearly with radius, while 
the negative chiral-cohesion energy, E tW i S t , both vanishes 
faster than linearly as R — > and diverges slower than 
linearly as R — > oo. Hence, the surface term generically 
dominates the bundle free energy in the small and large 
bundle limits. If the surface energy is above a critical 
value, S c , then the negative cohesive energy gain from 
bundle twist is insufficient to lower the free energy below 
the single-filament minimum at R = 0. In this case we 
have a state of dispersed single filaments in solution. It 
is straightforward to show that, 



S c (Ae = 0) 



3K 3 9f n , 
4A 3 j 



(22) 



6o(\j,\_) is the maximum twist angle of the 



where 9 n 
bundle according to eq. (fT9|) . 

'7t w A. 



/ 7twA 3 _l \ 
V K, ) 



1/3 



(23) 



When £ is decreased below this critical value, there is 
a global minimum at some Ro for which F(Rq) < 0, in- 
dicating a thermodynamically stable state of finite-sized 
bundles. Just at S = S c the equilibrium size remarkably 
depends only on the elastic resistance to bundle twist, 
Ro(T, c ) = X s ±. A simple estimate of this length scale 
may be deduced from the dimensional analysis of fila- 
ment arrays [l2| . K% is largely a measure of the intrinsic 
bending cost of bending individual filaments, suggesting 
K 3 ~ po^p, where £ p is the persistence length of the fil- 
aments. One the other hand the in-plane elastic moduli 
largely reflect forces between neighboring filaments, gov- 
erned largely microscopic distances of order do ~ p 1 ^ 2 . 
Hence, for a columnar array near the limit of thermal 
stability we expect — Pq 3 ^- These estimates sug- 
gest that A 3 j^ ~ (ipdo) 1 / 2 - As the persistence length 
of biological filaments may be on the order of microns, 
it is clear that the chiral assembly mechanism is quite 
consistent with the observation of mesoscopic bundles f- 
actin which are 10s of filaments in diameter [3l|. Note 




M II = 

dispersed " S * 
filaments 

\ 

ft] * 

dispersed, 
twisted bundles 



bulk, filament. 

aggregates 



Ae = Ae 

FIG. 4: The diagram of state for the assembly of chiral 
filaments into straight bundles of twisted filaments. When 
the net free energy gain per unit aggregated filament length, 
Ae < 0, filaments are dispersed in solution in one of two 
states. Above a critical value of surface energy, E c , filaments 
do not aggregate, and for E < E c bundles of finite radius are 
stable. The blue curve shows the boundary between these 
states for columnar, liquid-crystalline bundles and the red 
curve shows the same boundary for a strongly cross-linked, 
solid bundle. 



also that these scaling arguments suggest an estimate of 
the critical surface tension in eq. (f2"2"| in terms of the 
persistence length: E c ~ dQ 2 (£ p /d ) 1 ^ 2 9f nax (in units of 
ksT). When the surface energy is further lowered below 
S c ,the preferred bundle size grows, ultimately diverging 
as R ~ ST 3 . 

For the case when Ae < 0, finite-sized bundles may 
still be stable below a critical surface energy, albeit 
with diminished size. In the limit that — Ae is very 
large, the bundle size is reduced, ultimately falling of 
as, Rq ~ (— Ae) -3 / 2 . This result is obtained by bal- 
ancing the positive energy penalty, — 7ri? 2 poAe, with 
the negative energy gained by twisting very small bun- 
dle. This size scaling also suggests that the critical 
value of the surface tension falls as Ae — > — oo, since 
2E C = i? p Ae - E twlst /Ro ~ (-Ae)" 1 / 2 . Therefore, 
for sufficiently low surface energy and below the point 
of bulk condensation, this model predicts that a state of 
dispersed bundles will always be stable. This thermody- 
namic behavior is summarized in Figure [4] 

Finally, we note the E = behavior of bundles as Ae 
approached zero from below. In this limit, the bundle 
size is much greater than A 3 _l , with the cohesive penalty, 
7ri? 2 p|Ae|, balanced against the chiral cohesion energy 
in the shear-dominate, large-bundle regime, which grows 
as — i? 2 / 3 . Hence, as the point of bulk condensation is 
approached along the line of vanishing surface energy, 
the equilibrium bundle diverges as i?o ~ |Ae|~ 3 / 4 . In this 
way, the point S = and Ae = is like a critical point, 
through which the state of the system may transition 
continuously from dispersed, finite-sized bundles to bulk 
aggregates of unlimited size. 
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B. Solid Filament Bundles (/^y / 0) 

We now briefly consider the effect of a non-zero resis- 
tance to shears that slide the filaments with respect to 
one another along their axes. Clearly, shear resistance is 
strong for filamentous protein bundles that are strongly 
cross-linked by binding proteins. But we see here that 
even for the case of a very weak resistance to shear, i.e. 
Mn/M-L ^ 1' the qualitatively different mechanical re- 
sponse to twist produced by Tiii becomes relevant as the 
size a growing bundle diverges near the point the point 
E = and Ae = 0. 

A non-zero resistance to sliding shears alters the ther- 
modynamics of bundle assembly from the /zy = case 
described above in two important ways. First, the pres- 
ence of another positive elastic modulus in the model 
necessarily raises the free energy of a bundle configura- 
tion, with a given value of R and 9. Hence, a non-zero 
value of jiu , reduces the range of thermodynamic stability 
of finite sized bundles. The most straightforward anal- 
ysis of this occurs when 6n± is non-zero but sufficiently 
smaller to consider the effect of out-of-plane shear per- 
turbatively. For /q| ^ 0, we consider a more general twist 
energy, 



Ei 



1,R) 



4 



-2 7 Tw0i?, 



(24) 

where recall from (fTo) that 9 2 ± oc ^y//ij_. If 9 max ^> 
9, then the resistance to sliding shear only corrects the 
predictions of the /iy = of Rq and 9o at O[(#y_i_/#o) 2 ]- 
Therefore, to leading order the free energy correction per 
unit length from sliding shear is simply (K3/ 4)ir9 2 ± 9 2 l R 2 l . 
From the condition that F(Rq) — at £ c we can estimate 
the leading order reduction of the critical value of surface 
energy needed to stabilize dispersed bundles, 



AS,. 



K 3 Ro t 

8A 3_L 



9 2 



(25) 



where 9q and Ro are determined according the /zy = 
model. At the point of bulk filament condensation the 
relative size of this correction (from eq. (f2"2"|) is, 



/ A£ c 
V E r 



Ae=0 



69 2 ' 



(26) 



The predicted reduction of the thermodynamic stability 
of dispersed bundles in the presence of non-zero modulus 
flu in the diagram of state is shown in Figure [4j 

The second principle consequence of a resistance to 
sliding shears is the reduction in equilibrium size of dis- 
persed bundles. We consider this effect along the line 
Ae. As described above for the columnar the 
bundle surface energy decreases below £ c , the size of the 
bundle initially decreases as i?o ~ E_ 3 . According to 
eq. (|19p this is accompanied by a consequent drop in 

— 1/3 

bundle twist, 6q ~ R () . Therefore, even in the limit 



En — A, 




lnE 



s=s cOiN = o) 



FIG. 5: A plot of equilibrium radius of straight bundles of 
twisted filaments along the line of bulk condensation, where, 
Ae = 0. For a bundle that is strictly liquid crystalline, 
jUy = 0, the growth of the bundle is shown by the black curve. 
The equilibrium curves spanning the range of weakly solid, 



10" 



to strongly solid, 0u ± — 10 , shown in color. The 



increased resistance to out-of-plane shear reduces the bundle 
size, hence, Ro monotonically decreases with 0^ ± . This figure 
highlights the 6\ x crossover for the divergence of bundle 
size from columnar response (Ro ~ S" 3 ) to solid response 
(TJo-E" 1 . 



of vanishingly weak shear response, as the critical point 
(Ae = £ = 0) is approached the equilibrium bundle nec- 
essarily reaches a point where sliding shear dominates the 
elastic response, indicated by the fact that #o falls below 
0ii _i_ . This condition is met for sufficiently large bundles 
such that, Rq > A3_L(6' maa; /6'||j_) 3 . In this low surface 
energy limit, the sliding shear dominates the mechani- 
cal response to twisted bundle growth. Hence, equilib- 
rium bundle twist decays with R more rapidly than the 
/iy = case, 9q ~ R- Given this level of twist, from eq. 
(|24")l. the dominant size-dependence of E tW i S t results from 
the in-plane mechanical resistance which grows as 1/R 2 . 
Balancing this energy against the surface energy of the 
bundle, we find that the growth of the bundle generically 
crosses over to Rq ~ in the E — > limit. This is a 
notably weaker divergence than the Rq ~ E -3 growth of 
the columnar- hexagonal bundle model. Figure ?? shows 
equilibrium dependence of a solid filament bundle radius 
on the surface energy for a range of mechanical behavior, 
from /iy <C fix to A 1 !! ^ I 1 !-- 



IV. WRITHING BUNDLES OF HELICAL 
FILAMENTS 

In this section, we explore a separate microscopic 
model of filament chirality that also leads to the self- 
limited growth of bundles. Many biological filaments 
adopt naturally-helical configuration in their ground- 
state. The prototypical example of such filaments is the 
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bacterial flagella. Due to structural confirmations of the 
protein subunits composing flagella, these filaments are 
even known to adopt both right- and left-handed helical 
confirmations [l6|. Here, our task is to demonstrate that 
this type of intrinsic chiral structure also frustrates the 
assembly of densely-packed filament bundles. 

As a microscopic model of intrinsically helical fila- 
ments, we consider filaments that are anisotropic in their 
cross-section and have a tendency to bend around an axis 
that itself rotates around the filament tangent at a fixed 
rate. Following the standard linear elastic description of 
rod deformations (25j we introduce an orthonormal co- 
ordinate frame: ei x e 2 = £3 = t. A filament with a 
preferred helical configuration can be described by the 
following elastic energy, 

Hrod = 7; J ds [C Kl (Kl - K ) 2 + C Ki n\ + C W (CJ - Wq) 2 ]. 

(27) 

Here, J ds denotes an integral over filament arclength 
and Ki — e, ■ <9 s e3 describes bending in the two distinct 
in-plane directions and w = (d s e\ x ei) • e 3 describes the 
rod twist, or rate at which anisotropy of the cross-section 
rotates around the tangent direction. C\ , C2 and C3 are 
the bend and torsional moduli that penalize deformations 
from the ideal, helical state of the filament K\ — kq, k,i — 
and lo = luq- 

When assembled into a hexagonally-coordinated bun- 
dle, geometrical constraints make it impossible for fila- 
ments to maintain their ideal configuration. To capture 
the elastic cost of packing non-ideal filaments we intro- 
duce a specific parameterization, 

ex = cos ^>n + sin iph, (28) 

= — s'mtph + cos -0b, (29) 

where tp is a function of filament arc-length, and n and 
b are the unit normal and bi-normal of the Frenet frame 
of the filament backbone [29J . In these coordinates, the 
elastic energy of the filament can be computed in terms 
of tp and the geometry of the filament backbone, 

W TO d = ^ J ds [C K (k cos tp - k ) 2 + C K n 2 sin 2 tp 

+ CU^' + t-loo) 2 ], (30) 

where K and r are the respective curvature and torsion 
of the backbone curve, and we consider the simplified 
case C Kl = C K2 = C K . For a given backbone geome- 
try, we may compute the induced twist of the filament 
by solving the Euler-Lagrange equations for tp, which are 
in the most general case, non-linear and inhomogeneous. 
For the purposes of the following analyses, we will con- 
sider the two limiting cases. In the of easy twist, when 
C K C u , the filaments lock into the preferred state of 
bend: tp — 0. In the limit opposite limit of easy bend, or 
when C K -C C u , the groundstate achieves ideal torsion: 

tp' =LJ -T. 




Q=0 O = -T O = -2T 

"over" -twisted untwisted "under" -twisted 

FIG. 6: Three helical bundles described by eqs. (pTj) and J32}. 
An isometric bundle-with a constant separation between all 
nearest neighbor filaments-can only be obtained if there is no 
net twist, SI + T = (i.e. the coordinate frame compensates 
for the natural rotation of the Frenet frame). 

In the presence of strong adhesive forces, helical fil- 
aments form bundles with a preference for a uniform 
center-to-center spacing between filament backbones. 
Thus, as in the case of straight filaments, in the plane 
perpendicular to the backbone, the unstressed state of 
inter-filament forces is hexagonally packed. The geome- 
try of the bundle can be described in a similar manner 
to the single-filament geometry. The central axis of the 
bundle is described by a central curve, Ro(S'), with the 
associated Frenet frame, T, N and B. Note that the up- 
per case vectors refer to the geometry of bundle, while 
the lower case vectors refer to the geometry of a single fil- 
ament. Given the central curve, Ro, the positions of the 
filaments in the bundle correspond to two coordinates, 
x and y, in a coordinate frame that rotates along the 
bundle trajectory: 

X = cos^S^N + sin(05)B, (31) 

and 

Y = — sin(ClS)N + cos(nS)B. (32) 

These coordinates allow for an additional twist of fila- 
ments, n, relative to the natural rotation of the Frenet 
frame around the bundle axis. Using these coordinates 
the tangents of individual filaments can be derived, 

j = (Q + TH + (1-/Cr-N)T 
^J(n + T) 2 r 2 + (I - K, r-N) 2 

Here, r = xX + yY is radial separation of filament from 
the central axis of a bundle and (p = (—yX+xY) / r, is the 
direction of filament tilt in the plane of hexagonal order. 
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The respective curvature and torsion of the bundle axis 
are denoted by K, and T. 

According to the non-linear contributions to the in- 
plane strain tensor in eq. |1| and cq. ([55]) when 
Q + T ^ 0, 7^ 0, indicating a failure to maintain 
a constant separation between neighboring filaments in 
the assembly. In the case when 57 = — T, filaments are 
collinear; all filaments share the common tangent of the 
bundle axis, t = T. This is condition for an isometric 
packing of filaments [2l|, in which the interfilament sep- 
aration is constant in the bundle, throughout the cross- 
section and along the length. We should be clear to point 
out that these structures as distinct from the isometric 
textures considered by Achard et al. to model helical rib- 
bon formation in bent-core liquid crystals [HJ which can- 
not support a uniform hexagonally geometry in the cross- 
section. In the following analysis we focus on this limit of 
isometric packing favored by strong adhesive forces be- 
tween neighboring filaments in the bundle. An isometric 
bundle is untwisted (i.e. t • (V x t) = 0), distinguishing 
this mechanism for chiral filament assembly from the one 
that discussed in Sec. IIII1 

In an isometric packing all filaments share the Frenet 
frame of the bundle axis, so that differences in filament 
geometry derive from differences in arc-length element of 
filaments, ds, that of the central bundle axis, dS: 



ds=(l-tCr- N)dS. 



(34) 



From this relation, it is straightforward to compute the 
relationship between filament geometry and geometry of 
the writhing bundle, 

(35, 



l-K. r-N 



1-/Cr-N 



These relations reveal the geometric frustration inherent 
to periodid systems with preferred curvature: it is not 
possible to assemble constant-curvature filaments or lay- 
ers with a uniform separation [32l ? ]. As bundles grow 
to larger and larger radii, the curvature and torsion of 
individual filaments necessarily diverges. It is this frus- 
tration that may lead to a thermodynamic limit to lateral 
size of a bundle of helical filaments. 



A. "Easy Twist" Model (C K > C w ) 

We first consider the case of easy twist, in which we as- 
sume filaments always bend towards the preferred axis, 
ei< This case arises naturally in the limit that the bend- 
ing modulus is considerably larger than the torsional 
modulus of the filaments. Using the results for the lo- 
cal variation of filament bend and torsion and setting 
if> = in eq. (|30"|) . we derive the following elastic cost of 
an isometric filament packing, [52] 



(IC 2 C K + T 2 C^ 



1-y/l- (KR)< 



{ICR) 2 



-C K {K- k ) 2 R 2 + C U {T- cjofR 2 



for ip = 0. (36) 



Here, po is number of filaments per unit area in the cross 
section of the bundle. The terms in the first line of eq. 
(|45| account for the singular bending of filaments at the 
edge of the bundle in limit that K.R — > 1. This singular 
dependence of bending energy on R leads to a diverging 
lateral stress, restraining the lateral growth to R < Kr 1 . 

To analyze the thermodynamics of aggregation we con- 
sider the formation of bundles in the presence of cohesive 
interactions leading to net negative free energy contribu- 
tion per bundled filament, — irpoR 2 LAe. Combining the 
cohesive energy with the elastic cost of packing helical 
filaments, expanded to 0[i? 4 ] we obtain the free energy 
for growing bundles in the easy twist model, 



F(!C,R) 
npoL 



^(JC- Ko ) 2 -Ae}R 2 + ^(ICR) 4 



for ip = 0. (37) 



As we are working in the limit that C^/C K — > 0, we have 
dropped the terms of 0[C U }. 

When Ae > 0, filaments aggregate. However, these 
aggregates remain finite in radius below a critical value 
of the cohesive free energy per unit length of bundled 
filament, 



Ae r 



C K 



(38) 



In the range, < Ae < Ae c , bundles of finite radii are 
thermodynamically preferred, while for Ae > Ae c fila- 
ments unwind form bulk aggregates. 

Minimizing over bundle radii, we find a relation be- 
tween bundle size, Rq, and curvature, 



R 2 {K) = 



[2Ae-C K (lC-Kq) 



0. 



Ae < %*(/C - «o) 



Using this to rewrite the free energy in terms of bundle 
curvature only we find for Ae > 0, 



irp L 



-2 







[C K (fC-K, ) 2 -2Ae] 2 



Ae > ^f{K - k ) 



Ae< ^(JC-k ) 2 
(40) 



Finally, optimizing F in terms of equlibrium bundle cur- 
vature, /Co, we find, 



/C = Kofi 



Ae 
Ae7 



(41) 



From this result we see directly that aggregation leads to 
the straightening of filaments in the limit that Ae — > Ae c . 
In turn, this result yields the following prediction for the 
equilibrium size of bundles, 



Rn — K n 



2Ae 



(Ae c - Ae) 3 



(42) 



11 




FIG. 7: Plots of equilibrium structure of bundles of helical filaments for easy twist model (red curves) and easy bend model 
(blue curves). The dependence of the curvature of the central axis of the bundle, /Co, on the cohesive energy gain per unit 
filament length is shown in (a). The dependence equilibrium size of bundles on cohesive energy is shown in (b), with the inset 
highlighting the singular growth at the onset of bundling: _Ro ~ (Ae) 1 / 2 . The ratio of bundle torsion to bundle curvature, 
To //Co, is shown in (c). This ratio corresponds to the "aspect ratio" of the helical bundle, helical pitch relative to helical radius 
of bundle axis. 



Here, we find that upon aggregation the structure of the 
bundle evolves continuously from Ae = to Ae = Ae c . 
As the cohesive strength is increased, the equilibrium 
value of K, decreases from the preferred value, kq, to ac- 
commodate a greater number of filaments subject to the 
constraints that JCqRq < 1. Note, in particular, the sin- 
gular growth at the two limits of finite-sized aggregation: 
Ro ~ (Ae) 1 / 2 for Ae -> 0+ and i? ~ (Ae c - Ae)" 3 / 2 for 
Ae — * Ae c . These results are summarized in Figure [7] (a) 
and (b). 

A further measure of the structural evolution of bundle 
can be obtained by retaining terms linear in C u in eq. 
([57]) and calculating the evolution of the preferred bundle 
torsion, 



% = T 



2(Ae c - Ae) 
2Ae c - Ae 



(43) 



This predcits that individuals are untwisted just as they 
are unbent by stronger cohesive forces. The ratio Tq/JCo 
determines the "aspect ratio" of the helical bundle, that 
is, the ratio of the pitch to the radius of the path fol- 
lowed by the central axis of the helical bundle. The re- 
sults for Tq and /Co show explicitly that from the onset of 
aggregation, when Tq/K-q = tq/kq, to the point of bulk 
aggregation, when To //Co = 2t /k Qi the helical shape of 
the growing bundle continuously stretches along the pitch 
axis (see Fig. 0(c)). 

As in the case of straight bundles of chiral filaments, 
we should also attribute a positive energy cost, E, to the 
surface of bundle. Together with the adhesive energy 
gain per filament, Ae, E determines the thermodynamic 
stability of finite-radius aggregates of helical filaments. 
We estimate a critical value of surface energy, E c , above 
which finite-size aggregates are not stable from the con- 
dition 27ri?oiE c + F(/Co, Ro) = 0. This calculation yields 



dispersed, 
helical filaments 




bulk, straight 
aggregates 



/ Ae r 

dispersed, 

helical bundles 



Ae 



FIG. 8: The predicted diagram of state for the assembly of he- 
lical filaments in solution. If the net cohesive free-energy gain 
per bundled filament is below a critical value both dispersed 
filament and dispersed filaments may be stable, depending on 
the size of the surface energy of the bundles. The bound- 
aries between these to states are shown for the easy twist (red 
curve) and easy bend (blue curve) models. 



the following dependence of E c on cohesive energy, 

(Ae/Ae c ) 3 / 2 



E c (Ae) = C K n po 



VI - (Ae/Ae c ) 



(44) 



Because E c diverges as Ae — ► Ae c , the easy twist model 
predicts that as the system approaches the point of bulk 
aggregation from the state of dispersed single filaments, 
it necessarily passes state of dispersed bundles of finite 
diameter. The predicted diagram of state is shown in 
Figure [5] 
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B. "Easy Bend" Model (C K < C w ) 

We now consider the opposite case of easy bending, 
in which filaments lock into the the preferred state of 
twist despite the constraints imposed by isometric fil- 
ament packing in the bundle. From eq. (|30p . this is 
accomplished when d s ip — loq — r. Using the result 
of eq. (|35p for the variation in filament torsion and 
setting T — ujq, yields the following filament rotation 
ip = —nr sin(0 + ujqS). 



E m 6 

■np L 



= C K 



l-y/l- (ICR) 2 - 



(ICR) 2 



ICR Jx{\K\R) 



(ICR) 2 



+ ^(JC~k ) 2 R 2 



for ip = ujq 



(45) 



Note that in comparison to the C u = , easy-twist model 
in eq. I|36p . this result only differs by the term in square 
brackets on the second line. This is a signature of a 
more rapid increase of the excess bending induced by the 
bundle geometry in this limit 

Including the cohesive energy gain per unit length 
of filament and expanding this result in the limit that 
ICR <C 1 as before we obtain, 



F()C,R) 
npoL 



%(/C-Ko) 2 -A e li? 2 + %(/C 4 + Ko |/C| 3 )i? 4 
A Jo 



for ip' = uj 



(46) 



Again, notice that the correction due to imperfect fila- 
ment packing in the easy bend grows as |/C| 3 rather than 
the weaker /C 4 response of the easy twist model. Proceed- 
ing as before we minimize the free energy over bundle 
radii to find, 



, 2Ae-C fi (/C-Ko) 2 
C k (K 4 +k K 3 ) 



2Ae >C K {K- n f 
2 Ac < C K (IC ~ k q ) 2 



(47) 



Ae 1 / 2 , with the same exponent as the easy- twist model. 
Again, as cohesive strength increases and the size of the 
bundle grows, the bundle curvature is diminished, so that 
in the opposite limit, near the point of bulk condensation, 
/Co — > linearly as (Ae c — Ae) — * 0. Thus, as before 
when Ae — > Ae c , the equilibrium size of bundles diverges, 
albeit with a smaller exponent than in the easy twist case: 
Rq ~ (Ae c — Ae) -1 . Unlike the easy twist case, however, 
bundle torsion maintains the preferred values, Tq = u>o, 
for condition all cohesive strengths. Therefore, the aspect 
ratio of the helical bundle (pitch/radius) as characterized 
by 7o//Co, divergences as point of bulk condensation is 
reached, indicating the conformation of the bundle to be 
highly extended. This behavior is summarized in Figure 
d(a)-(c). 

Finally, we compute on estimate for the maximum sur- 
face energy, E c , for which finite-sized bundles are thermo- 
dynamically preferred over single helical filaments. The 
predicted diagram of state is shown in Figure [5J In con- 
trast to the easy twist model, we find a maximum critical 
surface tension at the point of bulk condensation deter- 
mined only by the preferred structure and bending elastic 
of a single filament: 



S c (Ae c ) 



C k kqPq 
3V3 



(49) 



Interestingly, in either the easy twist and easy bend lim- 
its, the assembly properties of the helical, isometric bun- 
dles is sensitive only the elastic cost of deforming fila- 
ments from their preferred state of bend, and not to C w 
or ro. From Fig. however, we see that the resistance 
to filament twist has important consequences for equi- 
librium assembly of helical filaments. Tuning the ratio 
C U /C K from (easy twist) to oo (easy bend) leads to a 
drastic reduction in the range of thermodynamic stability 
of dispersed bundles of finite diameter. 



CONCLUSION 



and, 



F(IC,Ro) 
npoL 



r, [C K (X-«o) 2 -2Ae] 2 
C K (/C 4 +K /C a ) 



2Ae > C K (K - /so) 



2Ae < C K (K, - K ) 2 

As before there is an upper limit to cohesive energy at 
which finite-sized bundles will form that is determined 
only by the energy cost per unit length of unbending 
the helical filaments: Ae c = G\n 2 ) /2. Unlike the easy 
twist case, it is not possible solve for preferred bend an- 
alytically. Nevertheless, it is possible to analyze bundle 
properties near the respective points of bundle formation 
(Ae — > 0) and of bulk aggregation (Ae — * Ae c ). 

It is straightforward to show that in the limit of Ae — > 
0+, /Co — kq minimizes F(JC, Rq). Hence, the equilibrium 
bundle size grows continuously from Ae = as Rq ~ 



In this study, we have explored two geometrically dis- 
tinct mechanisms by which chiral structure frustrates 
the two-dimensional assembly of filaments. In the first 
mechanism, forces between chiral molecules favor a pre- 
ferred relative twist of the orientation of neighboring 
molecules. While such a twist would be incompat- 
ible with a bulk system of hexagonally-ordered fila- 
ments, a globally twisted structure is compatible to two- 
dimensional packing in bundles of finite diameter. The 
competition between the effect of chiral interactions and 
the various mechanical costs of distorting the hexago- 
nal filament packing give rise size-dependence twisting of 
filaments around the central axis of the bundle. The size- 
dependent twist of the bundle, in turn, leads to a ther- 
modynamically preferred, finite lateral size of the bundle, 
provided that adhesive forces between filaments are suf- 
ficiently weak. 
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In the second mechanism, we consider the assembly of 
filaments whose preferred state is one of helical structure. 
The complexation of helical filaments induces a natural 
writhe to backbone of a growing bundle. Unlike the case 
of where interfilament forces induce a relative twist of 
molecular orientation, this mechanism allows filaments 
to form without an elastic cost for distorting the hexag- 
onal packing in the plane perpendicular to the filament 
tangents. The isometric packings preserve the nearest- 
neighbor separation between all filaments in the bundle. 
As a consequence of the preferred writhing the geometry 
of the bundle, in order for a bundle to grow filaments 
must be distorted, bent and twisted, from their preferred 
geometry, ultimately leading to a singular elastic cost as 
ICR — > 1. This elastic cost of distorting filaments pro- 
vides a thermodynamic limitation to the size growing 
bundles provided the adhesive energy gain per bundle 
filament is less than the mechanical cost required to un- 
bend a filament, Ae c . 

In many ways, eukaryotic organisms make their liv- 
ing through the constant assembly and disassembly of 
filamentous molecules. Functions as diverse as cell di- 
vision and locomotion are accomplished by means of 
relatively small class of filamentous proteins: f-actin, 
microtubules and intermediate filaments 33]. In vitro 
studies of the assembly behavior of filamentous proteins 
have contributed greatly to the understanding of physi- 
cal mechanism by which cells regulate cytoskeletal assem- 
bly and function [34j . In particular, it has been widely 
noted that when forces between filaments in solution are 
sufficiently attractive, dense bundles of aligned filaments 
form, with a limited diameter [HI [H, HE EE S H, . 
This is most surprising in view of apparent absence of 
long range forces between filaments in solution. Accord- 
ing to the thermodynamic description of such an assem- 
bly process, the classical nucleation model, there can be 
no equilibrium limitation to the growth in the presence 
of a net free energy gain per aggregated filament. To ex- 
plain this apparent contradiction, a number of theoreti- 
cal mechanisms for the self-limited growth have emerged, 
focusing variously on the specialized nature of forces be- 
tween filaments condensed in the presence of multi-valent 
counterions [4l|, [42j as well as the elastic cost of de- 
fects forced into bundles by rapid quenching [43[ or by 
the toroidal topology of bundles form by long strands 
of DNA 44j. In conflict with assumptions of "electro- 
static" mechanisms for limited growth, finite-sized bun- 
dles are observed even when bundles are condensed in 
the absence of multi-valent ions, for example, by de- 
pletion forces [37], [HI or through the incorporation of 
specialized cross- linking proteins [3^, H(| ■ Additionally, 
finite-diameter bundles of nominally straight filaments 
are readily observed, demonstrating that finite-bundle 
growth cannot as general rule be attributed to the com- 
plex topology of the bundle (4f| . 

In view of this phenomenology, it is quite natural to 
consider the chiral frustration of filament assemblies as a 
more generic mechanism for limiting bundle growth. Bio- 



logical filaments inherit chiral structure from the handed- 
ness of the subunits from which they are built, generically 
imbuing filaments with some measure of helical, screw- 
like structure. In this study, we have demonstrated the 
viability of this mechanism to limit the growth of chiral 
filament aggregates by way of a very generic model. Ap- 
plying the conclusions of the generic continuum-elastic 
description of hexagonally-packed bundles to a system 
of particular interest-that is, for biological filaments- 
requires a detailed accounting of the microscopic physics 
underlying the interactions between densely-packed fil- 
aments as well as the mechanical forces needed to dis- 
tort them. Certainly, some information about the phe- 
nomenological costs of the various geometrical distortions 
described in Sec. [IT] can be inferred from single- molecule 
mechanical measurements, or at the least from dimen- 
sional arguments. In particular, recent mechanical mea- 
surements of the elastic cost of bending (or unbending) a 
sinfle salmonella flagcllum allows us to estimate the range 
of cohesive forces for which bundles should be thermody- 
namically stable. The bending stiffness of these flagella 
is measure to be C K = 3.5pN /1m 2 , while the preferred 
curvature is in the range of kq ~ 1/xm -1 . From eq. (|38p 
we may estimate that bundles of these flagella are stable 
if the cohesive free energy per unit length is less than 
Ae c w fc^T/nm. 

The case of straight bundles of twisted filaments is 
more directly relevant to the bundle formation of filamen- 
tous proteins, such as f-actin or collagen. Here, it is nec- 
essary to determine the values of 7tw, & measure of the 
local preference for filament twist. Computing or mea- 
suring the strength of chiral forces in liquid crystalline is 
a notoriously difficult affair [l[ , largely due to the effect of 
positional and rotational fluctuations which considerably 
weaken the mutual torques exerted by chiral molecules. 
The two-dimensional order of bundle geometry simplifies 
this problem, as the in-plane positional as well as the 
orientational fluctuations of the filaments are restrained. 
The most detailed theory of the strength of chiral forces 
for biological filaments has been developed by Kornyshev, 
Leiken and coworkers This model treats biolog- 

ical filaments, like DNA, cylindrical rods, along which a 
helical charge pattern is distributed. In principle, this 
model could be adapted to compute the explicit depen- 
dence of the preference for a locally-braiding geometry 
of hexagonally-packed filaments on charge and structure 
of the helical geometry. However, it is unlikely the elec- 
trostatic forces alone determine the value of 7tw, since 
bundled protein filaments are brought to very close sepa- 
rations over order their molecular diameters and smaller. 
At this close range, interactions between biological fila- 
ments are likely very complex, including contributions 
from steric and hydrophobic forces between very het- 
erogeneous molecules. Indeed, atomic force microscopy 
observations of pairs of intertwined actin filaments [481 ] 
suggest that preference for twist for molecules in close 
contact is quite strong, as filaments readily wind helical 
around another on length scales much shorter than the 



14 



molecular persistence length, exposing a state of large 
bending stress which must be compensated by a lower- 
ing of the inter-filmament potential. Clearly, the develop- 
ment of a predictive and realistic model for the strength 
of twist-inducing forces in biological filaments is an out- 
standing problem with significant implications for the 
self-assembly properties of biological filaments. 
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APPENDIX: GEOMETRIC DECOMPOSITION 
OF CHIRAL ELASTIC ENERGY 

Here we demonstrate that the chiral elastic terms 
allowed by symmetry can be decomposed into two 
rotationally-invarient measures of the bundle geometry, 
the twist, or Tw, and the writhe, Wr, of the bun- 
dle [27|, 50] . For a bundle with filaments fixed at the ends 
or as closed bundle these two quantities are topologically 
constrained through the theorem, 



Lk = Tw + Wr, 



(A.l) 



where Lk is the linking number which counts the num- 
ber of links between the curves traced out by the bundle 
backbone and filaments in the bundle. In the present 
problem, there are no topological constraints that fix Lk, 
instead Tw and Wr appear independently as rotationally- 
invariant, chiral-symmetry breaking terms. 

Consider a filament in the bundle described by the 
curve, r(s), and the ath neighboring filament in the 
hexagonal array that is described by the curve r(s) + 
r a (s). The twist measures the local rotation of the of r a 
relative to r(s). The integrated twist of these two curves 
is given by, 

1 f t • (r Q x d s r a ) 
= ~2tt J dS ^ ' (A - 2) 

where t = d s r is the filament tangent. Note that 
ds r a = t a — t where t Q is the tangent to the neighbor- 
ing filament, which for smoothly- varying filament orien- 
tations can be approximated by t a ~ t + (r Q • V)t. Using 
this we compute the average twist by summing over the 



six neighbors in the hexagonal array, 

(Tw) = -i- [dsJ2*- [*a x (r a -V)t], (A.3) 

where we have divided by 6 to account for the number of 
filament pairs per hexagonal plaquette. Using the result 
that for a hexagonal lattice ^2 a (f a )i(f a )j — 3<5jj (where 
i and j specify directions in the plane perpendicular to 
t) we find the result that, 



2 ( Tw > = tt- / dst-(Vxt) 

Z7T 



(A.4) 



which is precisely the nematic twist of the director field. 
To lowest order in the displacement of filaments perpen- 
dicular to z this operator is given by, 



t • (V x t) ~ z ■ (V x d z u±). 



(A.5) 



It is straightforward to identify the chiral-energy coupling 
to twist by expanding terms in eq. ([6]) to first order 
in gradients of Ui to show that the elastic preference 
to twist is given by the sum of the three chiral elastic 
parameters, 7t w = 7 + l' + l" ■ 

The writhe of a curve is strictly a global measure of 
geometry [29(, requiring a description of the curve at 
distant points along the backbone. Using a theorem by 
Fuller [491, the writhe may be written as a local quan- 
tity by defining it relative to a reference curve. In this 
case, it is natural to choose the writhe-free z axis as a 
reference curve, as this represents the orientation of the 
bundle filaments in the undeformed reference state. In 
this case, the writhe of a filament may be written as an 
integral over a local quantity as, 



Wr 



1 

2^ 



ds- 



[t x (t • V)t] 
1 + t -z 



(A.6) 



This formula holds provided that everywhere 
along the contour of the curve. Again, by expanding the 
local writhe operator to lowest order in Uj_, 

5 • (t x d s t) 1„ ,„ „, . , „. 

1 + 1 • z £ 

so that we may associate the higher order contribution 
from chiral elastic term t • [(t • V) V x uj_] [from the third 
term in eq. ([5])] with filament writhe. From this expan- 
sion we deduce that the elastic preference for writhe is 
simply, 7wr = — l" ■ Note that relative to the chiral pref- 
erence for twist, the elastic coupling to writhe represents 
a higher-order term in uj^ and derivatives of u^, which 
we may regard as symptom of the "non-local" nature of 
writhe. 
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